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ABSTRACT 


A  high-resolution,  multi-level,  primitive  equation  ocean  model  is  used  to 
investigate  the  combined  role  of  seasonal  wind  forcing,  thermohaline  gradients,  and 
coastline  irregularities  on  the  formation  of  currents,  meanders,  eddies,  and  filaments  in 
the  California  Current  System  from  22.5°  N  to  47.5°  N.  An  investigation  of  the 
dynamical  reasons  for  the  generation  and  growth  of  meanders  and  eddies  is  conducted 
along  with  a  sensitivity  study  to  investigate  the  formation  of  the  Davidson  Current. 

Model  results  are  consistent  with  the  generation  of  eddies  from  instabilities  of  the 
equatorward  current  and  poleward  undercurrent  via  barotropic  and  baroclinic  instability 
processes.  The  meandering  equatorward  jet  south  of  Cape  Blanco  is  shown  to  be  a 
continuous  feature,  which  divides  coastally-influenced  water  from  water  of  offshore 
origin.  The  area  off  southern  Baja  is  shown  to  be  a  highly  dynamic  environment  for 
meanders,  filaments,  and  eddies,  while  the  area  off  Point  Eugenia  is  shown  to  be  a 
persistent  cyclonic  eddy  generation  region.  Both  the  Southern  California  Countercurrent 
rounding  Point  Conception  and  the  shoaling  of  the  poleward  undercurrent  are  shown  to 
play  important  roles  in  generating  the  Davidson  Current  in  the  fall. 
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I.  INTRODUCTION 


The  California  Current  System  (CCS)  is  the  classical  eastern  boundary  current 
(EEC)  system  located  off  the  west  coast  of  North  America  extending  approximately  from 
the  Strait  of  Juan  de  Fuca  southward  to  the  tip  of  the  Baja  Peninsula  (Hickey,  1998).  As 
in  other  EEC  regimes,  satellite  infrared  sea  surface  temperature  (SST)  imagery  of  the 
CCS  has  revealed  a  complex  flow  structure  with  seasonal  variations.  Dominant 
mechanisms  responsible  for  the  observed  large-scale  structure  within  the  CCS  have  been 
shown  to  be  seasonal  variations  in  alongshore  wind  stress  (Bakun  and  Nelson,  1991), 
coastline  irregularities  (e.g.,  Batteen,  1997),  bottom  topography  (e.g.,  Bceda  et  al,  1984), 
and  temperature  and  salinity  variations  (Batteen  et  al.,  1995;  Batteen  and  Vance,  1998). 

The  CCS  is  influenced  by  at  least  three  different  water  masses.  Cold,  low  salinity 
Pacific  Subarctic  Water  is  advected  from  the  north  while  warmer,  more  saline  North 
Pacific  Central  Water  enters  from  the  west.  Warm,  high  salinity  Southern  Water  is 
advected  from  the  south  by  the  poleward  undercurrent  (Hickey,  1998).  Unlike  other 
eastern  boundary  currents,  salinity  increases  rather  than  decreases  with  depth  in  the  CCS 
(Wooster  and  Reid,  1963). 

The  climatological  mean  CCS  is  comprised  of  several  large-scale  currents  (Figure 
1).  The  predominant  flow  is  the  California  Current  (CC),  which,  in  the  mean,  is  a  broad 

I 

(~10(X)  km  offshore),  shallow  (surface  to  -500  km),  relatively  slow  (-10  cm/s),  year- 
round  equatorward  current.  Embedded  within  this  basic  flow  are  the  California 
Undercurrent  (CUC),  the  Davidson  Current  (DC),  the  Southern  California  Countercurrent 
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(SCC),  and  the  Southern  California  Eddy  (SCE).  The  CUC  is  a  relatively  narrow  (~  10-40 
km  width),  weak  (-2-10  cm/s),  subsurface  poleward  flow  which  is  strongest  at  -100-300 
m  depth  and  varies  seasonally.  The  DC  is  a  weak,  inshore,  poleward  surface  flow  north 
of  Point  Conception  during  fall  and  winter.  Other  poleward  surface  flows  include  the 
SCC  south  of  Point  Conception  and  the  SCE  shoreward  of  the  Channel  Islands  within  the 
Southern  California  Bight  (Hickey,  1979, 1998). 

Recent  observational  studies  have  shown  that  the  CCS  is  not  the  stable  system  of 
currents  with  the  well-defined,  unchanging  stracture  suggested  by  Figure  1.  Rather, 
within  the  mean,  large-scale  structure  of  the  CCS,  there  exist  perturbations  in  the  form  of 
mesoscale  meanders,  eddies,  and  filaments,  which  vary  both  spatially  and  temporally 
(Bernstein  et  al.,  1977;  Chelton,  1984;  Stiub  et  al.,  1991).  Numerical  modeling  efforts  by 
Batteen  (1997)  simulated  this  complex  eddy-mean  flow  interaction  within  the  CCS  using 
a  high-resolution,  multi-level,  primitive  equation  (PE)  ocean  model  emphasizing  the  role 
of  wind-forcing  and  coastline  irregularities.  Specifically,  the  seasonal  variation  of  the 
alongshore  component  of  wind  stress  was  shown  to  be  critical  in  generating  realistic 
horizontal  and  vertical  structures  for  the  surface  equatorward  CC  and  the  subsurface 
poleward  CUC.  Both  baroclinic  and  barotropic  instability  mechanisms  contributed  to  the 
generation  of  meanders,  filaments,  and  eddies.  Additionally,  coastline  irregularities  were 
shown  to  help  “anchor”  upwelling  filaments  and  to  enhance  the  growth  of  meanders  and 
eddies.  In  particular,  the  region  off  Cape  Blanco  was  identified  as  the  location  where  the 
inshore  edge  of  the  CC  leaves  the  coast  and  develops  a  meandering  jet  to  the  south.  In 
recent  observational  (Batteen  et  al,  1995)  and  modeling  studies  (Batteen  and  Vance, 
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1998)  on  the  combined  effect  of  temperature  and  salinity  on  density,  it  was  shown  that  the 
distribution  of  salinity  as  well  as  temperature  is  important  in  defining  the  large-scale 
structure  and  circulation  of  the  CCS.  Although  the  modeling  study  of  Batteen  and  Vance 
(1998),  which  included  wind  forcing,  coastline  irregularities,  and  thermohaline  gradients, 
was  able  to  reproduce  some  of  the  currents  in  the  CCS  (i.e.,  the  CC  and  CUC),  the  model 
did  not  generate  the  DC,  SCC,  or  SCE  which  are  thought  to  originate  south  of  the  model 
domain  used  by  Batteen  and  Vance  (1998). 

In  this  study,  the  model  of  Batteen  (1997)  and  Batteen  and  Vance  (1998), 
originally  restricted  to  poleward  of  35°  N,  is  expanded  by  extending  the  model  domain 
southward  to  22.5°  N.  As  a  result,  the  entire  CCS  region  from  Baja  to  the  Washington- 
Canada  border  is  included,  which  allows  all  of  the  major  currents  of  the  CCS  shown  in 
Figure  1  as  well  as  filaments,  meanders,  and  eddies  to  be  generated.  Besides 
investigating  the  combined  role  of  seasonal  wind  forcing,  thermohaline  gradients,  and  a 
realistic  coastline  in  the  generation  of  these  features,  the  dynamical  reasons  for  the 
generation  and  growth  of  meanders  and  eddies  will  be  explored  using  the  energy  analysis 
technique  of  Batteen  et  al.  (1992). 

With  the  larger  model  domain,  two  key  issues  that  have  not  yet  been  resolved  can 
be  addressed  using  the  results  of  the  model  simulations.  The  first  issue  to  address  is  the 
seasonal  variability  off  the  Baja  peninsula,  which,  with  the  possible  exception  of  northern 
Baja  (which  has  been  regularly  sampled  by^CalCOFI  and  covers  the  region  north  of  Point 
Baja  to  the  U.S.-Mexico  border),  remains  a  data  sparse  region  with  a  poorly  known 
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seasonal  cycle  (Hickey,  1998).  The  second  issue  to  address  is  the  seasonal  appearance 
and  origin  of  the  DC,  which  is  not  fully  understood.  For  this,  in  a  sensitivity  study, 
results  from  a  second  experiment  with  a  straight  coastline  south  of  ~35®  N  will  be 
compared  to  the  basic  experiment  and  examined  to  determine  the  importance  of  coastline 
geometry  within  the  California  Bight  to  the  development  of  the  DC. 

This  study  is  organized  as  follows:  section  2  describes  the  numerical  model,  the 
forcing,  and  the  experimental  conditions  used  in  the  basic  experiment;  section  3  includes 
analyses  of  the  results  of  the  model  simulations  and  of  the  generation  mechanisms  for  the 
meanders  and  eddies  in  the  CCS;  section  4  utilizes  model  results  to  describe  the  seasonal 
variability  off  the  Baja  peninsula  and  the  seasonal  appearance  and  origin  of  the  DC. 
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II.  MODEL  DESCRIPTION 


A.  MODEL  EQUATIONS  AND  METHOD  OF  SOLUTION 


The  PE  model  in  this  study  was  originally  a  coarse  resolution  model  used  in 
closed  basin  studies  by  Haney  (1974).  It  has  recently'been  adapted  by  Batteen  (1997)  for 
eddy-resolving,  limited  EBC  regions  with  open  boundaries  to  the  north,  south,  and  west. 
The  model  is  multi-level,  non-adiabatic,  and  uses  the  beta-plane  approximation.  It  has 
both  baroclinic  and  barotropic  velocity  components  and  uses  the  hydrostatic  and 
Boussinesq  approximations  as  well  as  a  rigid  lid.  The  governing  equations  are  as 
follows: 


du  -I 
dt  Po  ^ 


UU  —  t  UU  .  ^  _4  „ 

~~r~~  T  2 


d^u 


ou  ov  ow  ^ 


-^=-pi 

01 


p=Po[i-tt(r-ro)+yfl(^-5o)] 


^=-A„v"r+/s:H-^ 
dt  ^ 


ds  d  S 

—  =  -A„W^S  +  K„  — 
dt  “  ^  oz^ 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 
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In  the  above  equations,  t  is  time  and  {x,y,:^  is  a  right-handed  Cartesian  coordinate  system 
with  X  pointing  toward  shore,  y  alongshore,  and  z  upward.  The  corresponding  velocity 
components  are  {u,v,w).  T,  S,  and  p  are  temperature,  salinity  and  pressure,  respectively, 
and  p  is  density.  Table  1  provides  a  list  of  other  symbols  found  in  the  model  equations, 
as  well  as  values  of  constants  used  throughout  the  study. 

A  space-staggered  B-scheme  (Arakawa  and  Lamb,  1977)  is  used  for  the  horizontal 
finite  differencing.  This  scheme  has  been  shown  by  Batteen  and  Han  (1981)  to  be 
appropriate  when  the  grid  spacing  is  approximately  the  same  order  as,  or  less  than,  the 
Rossby  radius  of  deformation.  The  horizontal  grid  spacing  is  9  km  in  the  east-west 
direction  and  11  km  in  the  north-south  direction,  while  the  internal  Rossby  radius  of 
deformation  is  ~30  km. 

The  model  uses  ten  vertical  layers,  with  constant  z-levels,  at  depths  of  13, 46,  98, 
182,  316,  529,  870,  1416,  2283,  and  3656  m.  This  vertical  scheme  is  designed  to 
concentrate  more  layers  above  the  thermocline  in  the  dynamically  active  portion  of  the 
ocean. 

The  model  domain  (Figure  1)  encompasses  the  west  coast  of  the  United  States  and 
the  Baja  Peninsula,  from  22.5°  N  to  47.5°  N  (2816  km  alongshore),  and  from  107.5°  W  to 
132.5°  W  offshore  (2304  km  cross-shore).  The  coeistal  boundary  of  the  model  domain  is 
closed,  and  has  both  tangential  and  normal  velocity  components  set  to  zero.  To  isolate 
the  role  of  wind  forcing  from  the  possible  coupled  effect  of  wind  forcing  and  topography 
(which  will  be  considered  in  a  separate  study),  bottom  topography  has  been  omitted  and 
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the  eastern  boundary  is  modeled  as  a  vertical  wall.  The  constant  depth  used  in  the  model 
is  4500  m.  A  modified  version  of  the  radiation  boundary  conditions  of  Camerlengo  and 
OBrien  (1980)  is  used  for  the  open  boundaries  to  the  north,  south,  and  west  with  some 
spatial  smoothing  applied  in  the  vicinity  of  these  boundaries  (Batteen,  1997). 

Biharmonic  lateral  heat  and  momentum  diffusion  is  used  in  the  model  with  the 
same  choice  of  coefficients  (i.e.,  2.0  xlO*^  cm^/s)  as  in  Batteen  (1997).  Holland  (1978) 
showed  that  biharmonic  diffusion  acts  predominantly  on  submesoscales,  while  Holland 
and  Batteen  (1986)  found  that  baroclinic  mesoscale  processes  can  be  damped  by 
Laplacian  lateral  heat  diffusion.  As  a  result,  the  use  of  biharmonic  lateral  diffusion 
should  allow  mesoscale  eddy  generation  via  barotropic  (horizontal  shear)  and/or 
baroclinic  (vertical  shear)  instability  mechanisms.  As  in  Batteen  (1997),  weak  (0.5  cm  /s) 
vertical  eddy  viscosities  and  conductivities  are  used  and  bottom  stress  is  parameterized  by 
a  simplified  quadratic  drag  law  (Weatherly,  1972). 

The  method  of  solution  is  straightforward  with  the  rigid  lid  and  flat  bottom 
assumptions  because  the  vertically  integrated  horizontal  velocity  is  subsequently 
nondivergent.  The  vertical  mean  flow  can  be  described  by  a  streamfunction  which  can  be 
predicted  from  the  vorticity  equation,  while  the  vertical  shear  currents  can  be  predicted 
after  the  mean  vertical  flow  is  subtracted  from  the  original  equations.  The  other  variables, 
i.e.,  temperature,  salinity,  vertical  velocity,  and  pressure,  can  be  explicitly  obtained  from 
the  thermod3ntiamic  energy  equation  (6),  salinity  equation  (7),  continuity  equation  (3),  and 
hydrostatic  equation  (4)  respectively  (for  more  complete  details  on  the  method  of  solution 
see  Batteen,  1997). 
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B.  FORCING  CONDITIONS  AND  EXPERIMENTAL  DESIGN 


Previous  experiments  by  Batteen  et  al  (1989)  investigated  the  role  of  steady 
alongshore,  upwelling-favorable  winds,  with  and  without  alongshore  variability,  on  the 
generation  of  eddies,  currents,  meanders,  and  filaments  in  the  CCS.  More  recently, 
Batteen  (1997)  studied  the  contributions  of  seasonal  wind  forcing  and  irregular  coastline 
geometry  to  the  generation  of  these  features  in  the  CCS  between  35°  N  and  47.5°  N. 
Following  Batteen  (1997),  seasonal  wind  forcing  and  irregular  coastline  geometry,  along 
with  thermohaline  gradients  as  used  by  Batteen  and  Vance  (1998),  will  be  used  to 
investigate  the  generation  of  similar  features  throughout  a  larger  domain  encompassing 
the  entire  CCS  (i.e.,  between  22.5°  N  and  47.5°  N). 

To  explore  the  effects  of  seasonal  wind  forcing,  the  model  is  forced  from  rest  with 
climatological  wind  fields  from  a  2.5  by  2.5  grid  of  European  Centre  for  Medium-Range 
Weather  Forecasts  (ECMWF)  near-surface  wind  analyses  (Trenberth  et  al,  1990).  The 
monthly  mean  stresses  based  on  twice  daily  wind  analyses  from  1980-1989  have  been 
interpolated  spatially  to  the  9  by  1 1  km  model  resolution  and  temporally  to  daily  wind 
values. 

Sample  wind  fields  used  in  the  basic  study  are  shown  in  Figure  2,  which  depicts 
the  seasonal  influence  and  migration  of  flow  around  the  North  Pacific  Subtropical  High 
over  the  area  encompassed  by  the  model  domain.  Figure  2a  shows  a  divergence  of  flow 
at  the  coast  in  the  vicinity  of  -40°  N  during  the  winter  as  the  winds  circulate 
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anticyclonically  around  the  Subtropical  High  to  the  south  and  cyclonically  around  the 
Aleutian  Low  to  the  north.  By  spring  (Figure  2b),  the  Subtropical  High  has  begun  its 
westward  and  northward  expansion,  producing  increased  equatorward  winds  over  most  of 
the  model  domain  and  causing  the  divergence  of  flow  to  move  north  to  ~44°  N,  off  the 
Washington  coast.  In  July  (Figure  2c),  equatorward,  upwelling-favorable  winds  dominate 
along  the  entire  coastline  as  the  Subtropical  High  reaches  its  maximum  extent.  Figure  2d 
shows  a  decrease  in  magnitude  of  equatorward  winds  during  fall  as  the  Subtropical  High 
begins  to  recede  southward. 

The  effects  of  thermohaline  gradients  in  the  CCS  are  included  using  monthly 
temperature  and  salinity  climatology  from  Levitus  et  al.  (1994)  and  Levitus  and  Boyer 
(1994)  to  initialize  the  model,  and  once  a  day  to  force  the  model  at  the  western  boundary 
(132.5°  W).  Forcing  occurs  over  the  upper  seven  levels,  which  are  initially  assumed  to  be 
zonally  homogeneous.  Since  the  lower  three  levels  do  not  exhibit  significant  horizontal 
variation  in  temperature  and  salinity,  constants  are  assumed  for  each  level.  Temperature 
values  used  for  levels  8  (1416  m),  9  (2283  m),  and  10  (3656  m)  are  2.56  C,  2.08  C,  and 
2.00  C,  respectively,  while  the  salinity  constant  used  for  the  lower  three  levels  is  34.7. 
Only  the  upper  layer  (13  m  and  46  m)  temperature  forcing  conditions  for  the  northern 
(47.5°  N)  and  southern  (22.5°  N)  boundaries  show  significant  seasonal  variability  with  a 
temperature  maximum  in  late  summer  to  early  fall  and  minimum  in  late  winter  to  early 
spring  throughout  the  region  (e.g..  Figure  3).  Below  these  depths,  both  the  seasonal 
temperature  fluctuations  and  temperature  gradient  weaken.  Conversely,  salinity  forcing 
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conditions  at  all  seven  levels,  which  show  less  (more)  saline  water  to  the  north  (south), 
have  no  significant  seasonal  cycle  (not  shown). 

The  design  of  the  study  is  as  follows.  To  focus  on  the  combined  role  of  seasonal 
wind  forcing  and  thermohaline  gradients,  model  integrations  start  from  a  state  of  rest  and 
once  a  day  the  model  is  updated  with  ECMWF  winds  and,  at  the  western  boundary  with 
temperature  and  salinity  fields.  In  the  basic  experiment,  a  realistic  North  American 
coastline  is  used  to  include  the  effects  of  irregular  coastline  geometry  over  the  entire  CCS 
region.  In  the  second  experiment,  a  straight  coastline  south  of  -35°  N  is  used  to  isolate 
the  effects  of  coastline  geometry  within  the  California  Bight.  To  allow  the  model  to  reach 
a  state  of  quasi-equilibrium,  the  basic  experiment  is  run  for  three  years. 

C.  ENERGY  ANALYSIS  TECHNIQUE 

The  energy  analysis  technique  used  is  the  same  as  that  used  and  described  by 
Batteen  et  al  (1992)  and  is  based  on  that  of  Han  (1975)  and  Semtner  and  Mintz  (1977). 
This  analysis  is  done  to  gain  a  better  understanding  of  the  types  of  energy  transfer  during 
unstable  flow  in  the  CCS.  A  brief  sununary  of  the  method  follows. 

Kinetic  energy  is  calculated  for  the  horizontal  components.  After  a  quasi-steady 
state  is  reached  where  the  total  kinetic  energy  is  nearly  constant,  mean  kinetic  energy 
(MKE)  and  eddy  kinetic  energy  (EKE)  are  calculated  using  the  sum  of  squared  mean  and 
horizontal  fields,  respectively.  Next,  the  available  potential  energy  is  calculated  and  used 
to  determine  when  a  quasi-steady  state  is  reached  and  when  statistics  should  be  collected. 
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Then  both  mean  and  eddy  available  energies  are  computed.  The  baroclinic  and  barotropic 
energy  transfers,  defined  by  Batteen  et  al.  (1992),  are  used  to  argue  for  the  type  of 
instability  mechanism  (e.g.,  baroclinic,  barotropic,  or  mixed)  leading  to  the  initial  eddy 
generation. 
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III.  RESULTS  OF  THE  BASIC  EXPERIMENT 


A.  SPIN-UP  PHASE 

Due  to  the  combination  of  thermohaline  gradients  and  wind  forcing,  different 
oceanic  responses  are  expected  depending  on  the  season.  In  the  winter,  in  the  poleward 
end  of  the  model  domain,  the  large  meridional  high  to  low  atmospheric  pressure  gradient 
due  to  the  decrease  of  temperature  poleward  (Figure  3)  establishes  an  onshore  geostrophic 
flow,  while  the  poleward  wind  stress  results  in  an  onshore  Ekman  flow.  On  approaching 
the  eastern  boundary,  the  onshore  flow  turns  and  forms  a  poleward  coastal  current  (e.g.. 
Figure  4a).  In  the  equatorward  end  of  the  model  domain,  the  smaller  pressure  gradient 
(due  to  the  smaller  warm  to  cold  temperature  gradient)  and  the  equatorward  wind  stress 
results  in  weak  onshore  geostrophic  flow,  offshore  Ekman  flow,  and  a  coastal 
equatorward  smface  current  (e.g..  Figure  4b) 

During  the  upwelling  season  (-April  to  September),  the  combination  of  a 
weakened  pressure  gradient  and  increased  equatorward  winds  over  the  entire  model 
domain  lead  to  a  strengthening  of  equatorward  flow  all  along  the  coast  and  a  weakening 
(strengthening)  of  onshore  (offshore)  flow.  An  equatorward  coastal  jet  is  discernible  as  a 
relatively  narrow  coastal  flow  in  the  north,  which  broadens  in  scope  as  it  progresses 
southward  (e.g..  Figures  5a  and  5b).  Maximum  speeds  in  the  northern  half  of  the  domain 
are  -20  cm/s  whereas  current  speeds  associated  with  cyclonic  eddy  activity  in  the  south 
are  -30-50  cm/s.  There  is  also  an  inshore  poleward  flow  discernible  at  -300  m  depth 
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with  maximum  speeds  of  ~  15-20  cm/s  which  can  be  traced  northward  across  the  entire 
model  domain  (not  shown). 

A  cross-section  of  meridional  velocity  at  day  195  (Figures  6b  and  6c)  shows  the 
typical  vertical  structure  of  the  equatorward  coastal  jet  and  poleward  undercurrent  during 
the  upwelling  season  from  the  northern  portion  of  the  domain  through  the  southern 
portion  of  the  domain.  In  the  northern  portion  (Figure  6b)  the  equatorward  coastal  jet 
axis  is  within  ~60  km  of  the  coast  and  extends  to  ~400  m  depth  with  a  surface  velocity  of 
~15  cm/s.  A  relatively  weaker  poleward  undercurrent  with  a  core  velocity  of  ~5  cm/s  at 
~300  m  depth  is  also  seen.  The  offshore  extent  of  the  core  of  the  undercurrent  is  confined 
to  ~50  km  of  the  coast.  In  the  southern  portion  of  the  domain  (Figure  6c),  the  coastal  jet 
axis  extends  to  -250  m  depth  with  a  surface  velocity  of  -25  cm/s.  A  poleward 
undercurrent  with  a  surface  velocity  of  -5  cm/s  is  found  within  -100  km  of  the  coast. 

As  the  coastal  jet  and  undercurrent  become  fully  established  (from  day  -195  to 
day  -285)  the  currents  become  unstable  and  form  meanders  (e.g..  Figures  5a-5b)  as  well 
as  cold,  upwelling  filaments.  As  the  meanders  intensify,  eddies  are  formed  in  the  coastal, 
equatorward  region  of  the  domain  (Figure  5b),  due  to  the  vertical  and  horizontal  shear 
instabilities  between  the  equatorward  jet  and  the  poleward  undercurrent.  The  eddies  are 
predominantly  cold  core  and  cyclonic,  on  the  order  of  100  km  in  diameter,  and  regularly 
extend  to  -50-100  km  off  the  coast.  Throughout  the  upwelling  season,  which 
corresponds  to  the  period  of  maximum  equatorward  wind  stress  (Figure  2),  meanders, 
filaments,  and  eddies  continue  to  develop  equatorward  of  -44°  N  (not  shown). 
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In  the  fall,  the  Aleutian  Low  begins  to  deepen  and  the  Subtropical  High  migrates 
toward  its  wintertime  position  in  the  south.  With  this  transition,  poleward  wind  stress 
returns  to  the  northern  portion  of  the  model  domain  as  the  regional  pressure  gradient 
strengthens.  The  poleward  flow  near  the  coast  fluctuates  in  depth  and  intensity  from  the 
northern  to  the  southern  portion  of  the  model  domain  (e.g..  Figures  7a-7b).  For  example, 
off  Cape  San  Lazaro,  at  -25°  N,  a  strong  poleward  undercurrent  is  discernible  below 
-300  m  depth  with  speeds  up  to  -15  cm/s  (Figure  7a),  while  farther  north,  off  Cape 
Blanco,  -43°  N,  poleward  flow  is  found  throughout  the  water  colunm  at  the  coast  with  a 
core  velocity  of  -7.5  cm/s  centered  at  -320  m  depth  (Figure  7b). 

B.  ANALYSIS  OF  EDDY  GENERATION  MECHANISMS 

The  dynamical  reasons  for  the  generation  and  stability  of  eddies  in  the  CCS  are 
examined  using  the  energy  techniques  used  and  described  by  Batteen  et  al.  (1992).  From 
the  energy  transfer  analysis,  the  location  and  magnitude  of  baroclinic  and  barotropic 
transfers  can  be  found  and  examined  to  argue  for  the  type  of  instability  mechanism  (e.g., 
barotropic,  baroclinic,  or  mixed)  which  leads  to  the  initial  eddy  generation.  Analysis  of 
typical  cross-section  of  the  meridional  velocity  component  (e.g..  Figures  6  and  7)  show 
that  there  is  considerable  horizontal  and  vertical  shear  in  the  upper  layer  currents  in  the 
CCS.  Barotropic  instability  can  result  from  horizontal  shear  in  the  currents,  while 
baroclinic  instability  can  result  from  vertical  shear  in  the  currents.  As  a  result,  both  types 


15 


of  instability  (mixed)  can  be  present.  Energy  transfer  calculations,  which  consist  of 
barotropic  (mean  kinetic  energy  to  eddy  kinetic  energy)  and  baroclinic  (mean  potential 
energy  to  eddy  potential  energy  to  eddy  kinetic  energy)  components,  were  performed  over 
the  time  period  when  the  meanders  and  eddies  developed. 

Results  from  the  instability  analysis  of  baroclinic  and  barotropic  transfers  for  days 
90  to  120  (not  shown)  show  the  presence  of  both  barotropic  and  baroclinic  instability 
mechanisms  (i.e.,  mixed).  Although  mixed  instability  mechanisms  are  present,  baroclinic 
instability  is  the  dominant  generation  mechanism  for  cyclonic  meander  development 
which,  ultimately  leads  to  the  formation  of  cyclonic  eddies  as  seen  north  of  Point  Eugenia 
by  day  195  (Figure  5b).  This  analysis  is  consistent  with  Figure  6a,  which  shows  a 
relatively  stronger  vertical  than  horizontal  gradient  in  upper  layer  current  shear  at  the  time 
of  the  initial  instability. 

The  instability  analysis  of  barotropic  and  baroclinic  transfers  for  days  180-195  in 
the  northern  portion  of  the  domain  (Plates  la  and  lb)  show  that  there  are  generally 
coincident  maxima  of  barotropic  and  baroclinic  instability  along  most  of  the  coast  from 
-35°  N  -  41°  N.  Note  that  the  scales  used  in  Plates  la  and  lb  are  different  in  order  to 
highlight  areas  of  relative  maxima.  These  regions  of  baroclinic  and  barotropic  (i.e., 
mixed)  instability  ordinarily  correspond  to  areas  of  meander  and  eddy  development  that 
occur  along  the  coast  by  day  285  (e.g..  Figure  8a).  It  can  be  seen  that  the  magnitude  for 
each  area  of  baroclinic  instability  is  generally  greater  than  its  coincident  area  of 
barotropic  instability,  indicating  a  dominant  baroclinic  energy  source  and  a  weaker 
barotropic  energy  source.  This  is  consistent  with  vertical  cross-sections  in  this  region 
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(e.g.,  Figure  6b)  which  shows  that  although  there  is  both  vertical  and  horizontal  shear,  the 
stronger  vertical  gradient  dominates  the  horizontal  gradient  in  the  upper  layers. 

The  instability  analysis  of  barotropic  and  baroclinic  transfers  for  days  180-195  in 
the  southern  portion  of  the  domain  (Plates  Ic  and  Id)  also  shows  coincident  maxima  of 
barotropic  and  baroclinic  instability  along  the  coast  extending  throughout  the  domain. 
Again,  the  scales  are  different  and  baroclinic  instability  once  again  dominates.  These 
results  are,  as  expected,  consistent  with  vertical  cross-sections  in  the  region  (e.g..  Figure 
6c)  which  show  that  strong  vertical  gradients,  in  general,  dominate  horizontal  gradients  in 
the  upper  layers.  An  exception  occurs  near  Point  Eugenia  (-28°  N)  where  the  magnitudes 
of  barotropic  instability  are  comparable  to  those  of  baroclinic  instability  which  results  in  a 
mixed  instability.  The  relatively  strong  horizontal  shear  off  Point  Eugenia  at  -28°  N, 
1 16°  W,  compared  to  regions  north  and  south,  can  be  seen  in  Figure  5b. 

The  instability  analysis  of  barotropic  and  baroclinic  transfers  for  days  270-285, 
near  the  end  of  the  upwelling  season,  in  both  the  northern  and  southern  portion  of  the 
domain  (Plates  2a-2d)  show  that,  although  both  baroclinic  and  barotropic  instability  are 
still  present,  barotropic  instability  dominates.  The  increase  in  barotropic  instability  can 
be  qualitatively  seen  by  comparing  Figure  5  with  Figure  8,  which  illustrates  that  a  general 
increase  in  horizontal  shear  from  days  195  to  285  has  occurred  throughout  the  domain.  In 
conclusion,  as  the  model  spins-up,  although  both  baroclinic  and  barotropic  instability 
mechanisms  are  present,  in  the  early  generation  period  baroclinic  instability  dominates 
while  barotropic  instability  is  the  dominant  mechanism  in  the  later  period. 
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C.  QUASI-EQUILffiRIUM  PHASE 


Longer  run  times  (~3  years)  with  model  output  fields  averaged  every  3  days  for 
the  months  of  January,  April,  July,  and  October  illustrate  the  seasonal  variability  and 
complex  structure  of  the  CCS.  By  year  3,  the  model  CC  core  takes  the  form  of  a 
meandering  jet  embedded  with  numerous  eddies  and  upwelling  filaments  (Plates  3a-3d). 
Although  the  mean  direction  of  flow  remains  equatorward  near  the  surface,  the 
circulation  contains  large  zonal  components  with  relatively  intense  onshore  and  offshore 
transports. 

Model  results  show  that  in  winter  (Plate  3a),  flow  at  the  coast  meanders 
equatorward  and  advects  colder  water  near  the  coast  offshore  and  warmer  water  from  the 
west  onshore  throughout  the  model  domain.  Both  cyclonic  and  anticyclonic  eddies  exist 
on  either  side  of  the  equatorward  coastal  jet  and  closed,  cold-core  eddies  on  the  order  of 
100  km  in  diameter  are  present  off  San  Francisco  at  ~38°  N,  126®  W,  and  in  the 
southwest  comer  of  the  domain  at  -25°  N,  129.5°  W.  The  coldest  water  found  in  the 
model  is  located  adjacent  to  the  southern  Baja  coastline  coincident  with  persistent 
northwesterly  winds  due  to  the  wintertime  position  of  the  Subtropical  High  (Figure  2a). 
Average  current  speeds  of  the  equatorward  surface  jet  range  from  -10-50  cm/s,  consistent 
with  observations  of  the  CC  using  drifters  (Davis,  1985). 

In  spring  (Plate  3b),  meander  activity  is  more  pronounced  along  the  model 
coastline  and  more  closed  eddies  exist  offshore  within  the  core  of  the  equatorward  model 
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CC.  Eddies  present  in  January  have  propagated  westward  at  speeds  of  ~  1-3  km/day, 
consistent  with  Rossby  wave  propagation  speeds.  The  coldest  water  has  moved  north  up 
the  Baja  and  Southern  California  coastline  in  conjunction  with  the  seasonal  migration  of 
the  Subtropical  High  and  subsequent  expansion  of  upwelling-favorable  winds.  Similarly, 
in  the  model,  spring  marks  an  increase  in  magnitude  of  the  equatorward  surface  jet  from 
speeds  present  in  January.  Using  CalCOFI  data,  Lynn  and  Simpson  (1987)  noted  that  a 
large-scale  intensification  in  equatorward  flow  along  much  of  the  coastal  CCS 
accompanied  the  onset  of  springtime  conditions. 

In  summer  (Plate  3c),  temperatures  adjacent  to  the  coast  in  the  northern  portion  of 
the  model  (i.e.,  above  -40°  N)  have  decreased  as  upwelling-favorable  winds  now 
approximately  parallel  the  coastline  throughout  the  domain  (Figure  2c).  Evidence  of 
cold,  offshore-flowing,  upwelling  filaments  exists  in  the  vicinity  of  Cape  Blanco,  Cape 
Mendocino,  south  of  Point  Sur,  in  the  Southern  California  Bight,  and  below  Point 
Eugenia.  These  filaments  extend  -80-  200  km  offshore  and  merge  with  the  meandering 
equatorward  jet  that  has  speeds  of  -30-50  cm/s  and  alongshore  wavelengths  of  -100-300 
km,  consistent  with  Brink  and  Cowles  (1991).  An  example  of  a  mesoscale  disturbance 
embedded  within  the  CCS  is  illustrated  by  the  cyclonic  meander  that  forms  off  Point 
Eugenia  in  January  (Plate  3a),  intensifies  in  April  (Plate  3b),  and  breaks  away  by  July 
(Plate  3c)  to  propagate  southwestward  in  the  CC.  Bernstein  et  al  (1977)  noted  similar 
mesoscale  activity  off  Point  Eugenia  using  data  from  CalCOFI  surveys  of  the  CCS  taken 
from  April  to  July,  1952.  In  the  vicinity  of  the  U.S.-Mexico  border  (-32.5°  N)  within  the 
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Southern  California  Bight,  a  division  of  flow  is  evident,  consistent  with  Reid  (1963),  in 
the  equatorward  jet  as  a  portion  of  the  jet  turns  shoreward  (Plate  3c).  This  division  of 
flow  is  consistent  with  the  summertime  formation  of  the  SCE  (Lynn  and  Simpson,  1987). 

The  spring  and  summer  results  are  in  good  agreement  with  recent  satellite  imagery 
(Strab  et  ah,  1991;  Strub  and  James,  1995),  field  studies  using  Lagrangian  drifters  (Barth 
et  al.,  1994;  Barth  and  Smith,  1996a,b),  and  with  the  hypothesis  that  the  offshore 
separation  of  the  equatorward  coastal  jet  in  the  vicinity  of  Cape  Blanco  marks  the  start  of 
a  continuous  meander  along  a  temperature  front  that  can  be  traced  equatorward 
throughout  the  entire  CCS  during  spring  and  summer.  In  particular,  both  the  model 
results  and  Strub  and  James  (1995)  (hereafter  referred  to  as  SJ95)  show  the  connection 
between  the  equatorward  flow  off  Oregon  and  the  jet  that  meanders  along  the  temperature 
front  off  California,  i.e.,  the  equatorward  flow  off  Oregon  and  the  meandering  jet  are  part 
of  a  continuous  flow  that  originates  from  farther  north.  In  particular,  SJ95  (see  Figure  2a 
of  SJ95)  and  the  model  results  (Plate  3c)  both  show  a  narrow  jet  that  is  found  on  the 
offshore  side  of  the  temperature  front  at  ~20-50  km  off  central  Oregon  at  -44.5°  N. 
Additionally,  the  similar  presence  of  coastal  equatorward  flow  leaves  the  coast  in  the 
vicinity  of  Cape  Blanco.  Downstream  from  Cape  Blanco,  both  results  show  the 
equatorward  jet  and  temperature  front  that  form  a  large  meander  around  Cape  Mendocino 
and  then  flows  southward  toward  the  coast  where  it  turns  south-southwest  near  Point 
Arena.  Both  results  also  show  the  jet  and  its  associated  temperature  front  between  Cape 
Mendocino  and  Point  Arena  formed  inshore  of  -125°  W. 
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By  early  fall  (Plate  3d),  as  the  upwelling  system  begins  to  weaken  with  the 
equatorward  migration  of  the  Subtropical  High,  the  sharp  nearshore  temperature  front 
present  during  summer  north  of  Cape  Blanco  has  become  more  diffuse.  The  coastal 
equatorward  surface  jet  in  this  region  is  now  oriented  farther  offshore  and  embedded 
streams  of  colder  water  flow  southwestward  into  a  cold-core  eddy  west  of  Cape 
Mendocino.  South  of  Point  Conception,  the  division  of  flow  is  more  pronounced  as  a 
branch  of  the  equatorward  surface  jet  rotates  onshore  and  poleward  while  the  main 
portion  of  the  jet  continues  southward  offshore.  This  flow  structure  within  the  Southern 
California  Bight  coincides  with  a  seasonal  maximum  in  the  SCE  from  summer  to  early 
fall  (Hickey,  1998).  As  in  January  (Plate  3a),  in  the  vicinity  of  Point  Eugenia,  a  cyclonic 
meander  in  the  equatorward  jet  is  discernible.  In  time  (not  shown),  this  meander 
intensifies,  closes  off,  and  propagates  southwestward  through  the  stream. 


D.  ANALYSIS  OF  ENERGETICS 

Horizontal  maps  of  the  upper  layer  mean  kinetic  energy  (MKE)  and  eddy  kinetic 
energy  (EKE),  averaged  monthly  in  the  quasi-equilibrium  phase,  are  shown  in  Figures  9 
(11)  and  10  (12)  in  the  summer  and  early  fall  in  the  northern  (southern)  portion  of  the 
model  domain.  Maps  of  MKE  and  EKE  are  suggestive  of  where  the  mean  and  eddy 
energy  sources  are  to  be  found  (Holland  et  al.,  1983). 
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A  comparison  of  MKE  (Figures  9  and  11)  with  the  mean  velocity  field  (Plate  3c 
and  3d)  shows  that  high  values  of  MKE  are  found  along  the  axis  of  the  meandering 
equatorward  jet,  in  the  coastal  regions  where  the  coastal  current  is  the  strongest,  and  in 
regions  of  cyclonic  eddies  offshore  of  the  jet  at  -41°  N,  128°  W,  -38°  N,  129°  W,  -34° 
N,  124°  W,  and  -28°  N,  125°  W.  Maximum  values  of  MKE  are  found  in  regions  where 
mean  equatorward  flow  on  the  eastern  (western)  side  of  anticyclonic  (cyclonic)  eddies 
enhances  the  flow  between  cyclonic  and  anticyclonic  eddy  pairs.  An  example  of  this  can 
be  seen  by  comparing  the  eddy  pair  at  -40.5°  N  126°  W  off  Cape  Mendocino  (Plate  3c) 
with  the  corresponding  MKE  maximum  (Figure  9a).  An  additional  example  of  this  can 
be  seen  by  comparing  the  eddy  pair  at  -25°  N,  1 19.5°  W  off  Cape  San  Lazaro  (Plate  3c) 
with  the  corresponding  MKE  maximum  (Figure  11a).  MKE  maxima  can  also  be  seen  in 
regions  of  cyclonic  eddies  offshore  of  the  jet.  For  example.  Plates  3c  and  3d  show 
increased  current  velocities  on  the  western  sides  of  the  offshore  cyclonic  eddies  at  -41.5° 
N,  128.5°  W  and  at  -38°  N,  129°  W  which  correspond  to  areas  of  higher  MKE  in  Figures 
9a  and  9b. 

Within  the  northern  portion  of  the  domain,  a  comparison  of  EKE  (Figure  10)  with 
the  mean  velocity  field  (Plate  3c  and  3d)  shows  high  values  of  EKE  where  the  jet 
meanders  off  Cape  Blanco  at  -43°  N,  126°  W  and  on  the  eastern  and  western  sides  of 
offshore  cyclonic  eddies  at  -41°  N,  128°  W  and  at  -38°  N,  128.5°  W.  ffigh  EKE  values 
are  also  located  along  the  axis  of  the  coastal  equatorward  jet  at  -45°  N,  124.5°  W  and  in 
the  region  of  the  inshore  cyclonic  eddy  at  -38°  N,  124°  W.  Within  the  southern  domain, 
a  comparison  of  EKE  (Figure  12)  with  the  mean  velocity  field  (Plates  3c  and  3d)  shows 
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maximum  FKF.  values  at  33°  N,  118°  W  where  the  jet  meanders  toward  the  Southern 
California  Bight,  in  the  vicinity  of  -28.5°  N,  1 16°  W  where  a  cyclonic  eddy  is  generated 
off  Point  Baja  between  September  and  October,  and  farther  off  Point  Baja  at  -30.5°  N, 
120°  W  where  another  cyclonic  eddy  is  generated  between  September  and  October. 
Regions  of  high  model  EKE  values  correspond  to  areas  where  eddies  are  likely  to  be 
generated,  such  as  off  Cape  Blanco,  offshore  and  downstream  of  Cape  Mendocino,  Point 
Arena,  and  Point  Baja,  in  the  Southern  California  Bight,  and  in  the  coastal  indentations 
on  either  side  of  Point  Eugenia.  A  time  series  of  EKE  (Figures  10a- 10c),  depicts  higher 
values  of  EKE  present  during  late  summer.  For  example,  the  maximum  values  of  EKE 
for  the  eddy  located  at  -38°  N,  129°  W  in  July  increased  from  -25-50  cmVs^  (Figure 

10a),  to  -75-100  cmW  by  September  (Figure  10b). 

A  comparison  of  MKE  (Figures  9  and  11)  with  EKE  (Figures  10  and  12)  show 
that  maximum  MKE  and  EKE  values  generally  occur  in  the  same  vicinity.  Additionally, 
for  the  entire  model  domain,  the  magnitudes  of  MKE  are  larger  than  the  magnitudes  of 
EKE  for  all  months  (including  those  not  shown).  This  is  consistent  with  the  results  of  the 
model  simulations  which  showed  that  eddies  are  generated  from  instabilities  of  the  mean 
CC  and  CUC  via  baroclinic  and/or  barotropic  instability  processes. 

Kelly  et  al.  (1998)  recently  examined  altimeter,  moored  ADCP,  and  drifter  data 
for  the  CCS  to  determine  the  spatial  and  temporal  stmcture  of  eddy  fields  and  to 
investigate  the  variability  of  near-surface  EKE  in  the  CCS  west  of  124°  W  between  33° 
N-40.5°  N.  A  qualitative  and  quantitative  comparison  of  EKE  results  from  the  present 
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study  with  their  observations  shows  good  correlation.  Both  Kelly  et  al.  (1998)  and  the 
model  results  show  higher  values  of  EKE  in  late  summer  to  early  fall  between  -36°  N- 
40°  N  and  west  of  -124°  W,  coincident  with  increased  equatorward  flow,  along  with 
minimum  values  in  the  spring  (not  shown).  Specifically,  maximum  (minimum)  RKF. 
values  observed  by  Kelly  et  al.  (1998)  within  their  study  area  were  -0.02  mVs^  (0.01 
mVs^)  while  maximum  (minimum)  values  in  the  model  were  -0.03  mVs^  (0.01  mVs^). 
These  values  are  also  similar  to  Strub  et  al.  (1997).  Additionally,  both  studies  reveal  a 
tendency  of  EKE  maxima  to  propagate  westward  in  time.  For  example,  the  F.KF. 
maximum  located  at  -38°  N,  128°  W  in  July  propagates  to  -38°  N,  130°  W  by  October. 
Thus,  the  results  of  the  model  instability  analysis  are  consistent  with  Kelly  et  al.  (1998) 
who  also  suggests  that  the  baroclinic  instability  mechanism  is  the  dominant  forcing 
mechanism  in  the  generation  of  eddies. 
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IV.  MODEL  SIMULATION  APPLICATIONS 


An  additional  goal  of  this  study  is  to  utilize  model 


issues  that  have  not  yet  been  resolved.  The  first 
Baja  peninsula  which,  remains  a 


results  to  address  two  key 


issue  is  the  seasonal  variability  off  the 


sparse  region  with  a  poorly  known  seasonal  cycle 
(Hickey,  1998).  TT,e  secoed  issue  Co  address  is  flee  seasonal  appearance  and  origin  of  the 

DC  Which,  is  not  fully  understood.  For  the  second  issue,  the  tesults  of  a  sensitivity  shtdy 
will  be  used. 


SEASONAL  VARIABILITY  OFF  THE  BAJA  PENINSULA 


Lynn  and  Simpson  (1987)  have  described  the  flow  and  water  mass  chamcteristics 
off  the  Baja  Peninsula  based  on  analyses  of  23  years  of  CalcOH  data.  Surface 
q  rwaid  flow  is  strongest  along  the  coast,  showing  significant  seasonal  variability, 
with  maximum  flow  of  -20  cm/s  occurring  in  March  and  April.  This  inshore 
equatorward  flow  is  pmsent  most  of  the  year  with  no  spring  transition  as  obse^ed  north 
of  the  border.  The  near-surface  water  mass  characteristics  described  by  Lynn  and 
Simpson  (1987)  include  a  low-salinity  minimum  a,  -50-100  m  depth  (Reid  «  al..  1958) 
close  to  the  coast  which  indicates  the  core  of  the  CC. 

In  the  model,  the  flow  results  show,  consistent  with  Lynn  and  Simpson  (1987). 
te  in  spring  and  summer  (e.g.,  Plams  3b  and  3c),  the  equatorwaM  surface  jet  off  Baja  is 
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and  CO.  .0  d.0  coas,  while  during  tal,  and  win«  (e.g.,  Pla«s  3a  and  3d)  fte  jei 
weakens,  meanders  farmer  offshore,  and  is  embeds  in  a  more  chaoric  flow  regime.  Ttre 
model  also  supports  maximum  equarorward  current  speeds  of  -20-30  cm/s  offshore 
during  spring  (e.g.,  Plate  3b),  coincident  with  the  strongest  equatorwaid  wind  speeds 
(Hgure  2b).  The  water  mass  characteristics  of  the  model  simulation  show,  consistent 
with  Lynn  and  Simpson  (1987)  and  Reid  e,  al.  (1958).  that  south  of  Point  Eugenia  a  low 

salinity  minimum  of  -33.7  (Figure  13)  is  present  to  -50  m  depth. 

The  region  southeast  of  Point  Eugenia,  which  reptesents  the  largest  coastline 

pemirbation  along  the  Baja  Peninsula,  is  the  locmion  of  a  seasonally  intensifying  cyclonic 
eddy  such  that  return  flow  north  of  fhe  eddy  center  transport  warmer  water  hundreds  of 
kilometers  offshore  between  October  and  January  (Lynn  and  Simpson.  1987).  The  model 
shows  that  Point  Eugenia  and  the  mgion  to  the  south  are  areas  of  snong  seasonal 
variability  as  weU  as  a  generaUon  region  for  cyclonic  eddies.  In  January  of  year  3  (Plate 
3a).  an  inshore  branch  of  the  equatorward  surface  jet  wifl.  speeds  as  high  as  -25  cm/s 
turns  away  from  die  co^t  in  the  vicinity  of  Point  Eugenia  to  rejoin  the  main  stmam 
offshore,  indicating  development  of  a  cyclonic  eddy,  m  April  (Plate  3b).  die  cyclonic 
meander  is  better  developed  and  has  begun  to  detach  from  the  vicinity  of  Point  Eugenin 
By  July  (Plate  3c),  the  meander  has  closed  off  to  form  a  cold-core  eddy  which  has 
propagated  southwestward  and  Is  now  centered  at  -25”  N,  1 18.5  ”  W.  A  thin  filament  of 
colder  water  extends  south  from  Point  Eugenia,  transported  by  the  high-speed  inshore 
near-snrface  equatorward  flow.  Additionally,  cyclonic  turning  is  once  again  evident  in 
the  equatorward  surface  jet,  but  at  a  slighfly  different  locadon  from  January,  west- 
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northwest  of  Point  Eugenia.  By  October  (Plate  3d),  cyclonic  turning  to  the  west- 
northwest  has  further  intensified  and  approaches  the  pattern  shown  for  January  (Plate  3a). 
Note  that  in  January  (Plate  la),  a  cyclonic  eddy  centered  at  ~26°  N,  123  W  has  also 
formed  and  propagated  westward  in  time  from  the  same  generation  area.  Thus,  according 
to  model  results.  Point  Eugenia  appears  to  play  a  significant  role  in  cyclonic  eddy 
generation  and  subsequent  westward  propagation. 

Although  a  data  sparse  region,  the  area  off  southern  Baja  is  believed  to  be  a  highly 
dynamic  environment  for  meanders,  filaments,  and  eddies  (Lynn  and  Simpson,  1987; 
Poulain  and  Niiler,  1989;  Poulain,  1990;  Niiler  et  al,  1989).  Model  results  from  year  3 
support  this  description  for  the  region  south  of  Point  Eugenia.  Shoreward  anticyclonic 
meanders  in  the  equatorward  surface  jet  (Plates  3a  and  3b),  sharp  temperature  fronts 
alongshore  and  offshore-flowing  cold  filaments  (Plates  3c  and  3d),  combined  with  the 
flow  field  effects  of  cyclonic  eddies  detaching  from  Point  Eugenia  to  propagate 
southwestward,  all  contribute  toward  a  highly  variable  and  complicated  current  structure 
in  the  model. 

An  example  of  near-surface  and  subsurface  variability  in  the  model  for  this  area  is 
illustrated  through  cross-sections  of  mean  meridional  velocity  taken  at  25°  N,  in  the 
vicinity  of  Cape  San  Lazaro.  In  January  (Figure  14a),  equatorward  flow  persists  near  the 
surface  with  a  maximum  of  ~25  cm/s  nearshore  and  a  second  maximum  of  ~35  cm/s 
farther  offshore.  A  poleward  undercurrent  structure  is  present  within  ~60  km  of  the  coast 
with  a  tnaYimiim  of  -10  cm/s  below  -300  m  depth.  In  April  (Figure  14b),  equatorward 
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flow  is  still  present  nearshore  but  adjacent  poleward  flow  offshore  is  now  surface- 
intensified  to  ~25  cm/s  and  extends  throughout  the  upper  water  column,  indicating  the 
presence  of  an  anticyclonic  eddy  nearshore.  Farther  offshore,  the  equatorward  jet  has 
intensified  to  ~50  cm/s  in  response  to  the  occurrence  of  maximum  equatorward  winds  in 
spring.  By  July  (Figure  14c),  alternating  poleward  and  equatorward  flow  patterns  are 
discemable,  which  are  consistent  with  the  presence  of  a  cyclonic  eddy  near  the  coast  and 
an  anticyclonic  eddy  offshore.  By  October  (Figure  14d),  the  nearshore  cyclonic  rotation 
has  migrated  slightly  westward,  giving  way  to  a  narrow  band  of  surface-intensified 
equatorward  flow  of  ~  1 5  cm/s  alongshore. 

B.  SEASONAL  APPEARANCE  AND  ORIGIN  OF  THE  DAVIDSON  CURRENT 

Chelton  (1984)  used  23  years  of  hydrographic  measurements  off  Point  Sur  and 
Point  Conception  to  show  that  in  both  locations  the  average  poleward  geostrophic  flow 
extends  to  the  surface  from  -October  to  February.  Lynn  and  Simpson  (1987)  refer  to  the 
seasonal  change  in  direction  of  the  average  surface  geostrophic  flow,  from  equatorward  to 
poleward,  that  occurs  within  -150  km  of  the  coast  during  fall  and  winter  as  the  Inshore 
Countercurrent  (IC)  instead  of  the  DC  and  suggest  that  the  appearance  in  fall  and  winter 
of  the  IC  is  associated  with  the  shoaling  of  the  CUC.  According  to  Hickey  (1979,  1998), 
this  poleward  flow  can  be  continuous  around  Point  Conception.  She  notes  that  south  of 
Point  Conception,  a  branch  of  the  equatorward  CC  turns  shoreward  then  poleward  and 
calls  this  flow  the  SCC  if  it  rounds  Point  Conception  and  the  SCE  during  periods  when  it 
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exists  as  a  closed  circulation  within  the  Southern  California  Bight.  She  also  points  out 
that  the  SCC  has  a  seasonal  maximum  which  coincides  with  the  development  of  the  DC 
north  of  Point  Conception,  but  that  evidence  also  exists  which  attributes  the  DC  to  a 
shoaling  of  the  CUC  (Pavlova,  1966;  Huyer  and  Smith,  1974). 

The  results  of  the  basic  experiment  show  that,  consistent  with  Chelton  (1984),  off 
both  Point  Sur  (-36.5°  N)  and  Point  Conception  (-34.6°  N)  the  average  poleward  flow  at 
the  coast  extends  throughout  the  upper  water  column  from  -October  to  February  (e.g.. 
Figures  15a  and  15b).  The  results  also  show,  consistent  with  Lynn  and  Simpson  (1987), 
the  seasonal  transition  of  the  poleward  undercurrent,  from  weak  subsurface  flow 
nearshore,  to  strengthening  and  shoaling  surface  poleward  flow,  as  illustrated  at  32°  N  in 
Figures  16a-16c.  In  September  (Figure  16a),  weak  (-2.5  cm/s)  poleward  flow  is  centered 
at  -175  m  depth  and  underlies  equatorward  flow  at  the  surface.  By  October  (Figure  16b), 
poleward  flow  extends  to  the  surface  with  a  core  maximum  of  -5  cm/s  centered  at  -40  m 
depth.  By  November  (Figure  16c),  a  relatively  strong  coastal  surface  poleward  current  of 
-12.5  cm/s  is  discernible. 

There  is  also  appreciable  evidence  in  the  model  to  support  a  continuous  surface 
poleward  flow  (i.e.,  the  SCC  as  illustrated  in  Figure  1)  from  the  Southern  California 
Bight  region  to  north  of  Point  Conception,  which  according  to  Hickey  (1998)  is  caused  by 
a  wrap-around  of  the  divergent  branch  of  the  CC.  hi  particular,  inshore  poleward  flow  is 
shown  to  be  continuous  around  Point  Conception  as  illustrated  for  November  in  Figures 
17a  (south  of  Point  Conception  in  the  Southern  California  Bight  region),  17b  (off  Point 
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Conception),  and  17c  (north  of  Point  Conception).  Thus,  the  basic  experiment  was  able 
to  generate  surface  poleward  flow  (i.e.,  the  DC),  which  can  be  attributed  to  both  the 
shoaling  of  the  CUC  (as  in  Figure  16)  and  to  SCC  rounding  Point  Conception  (as  in 
Figure  17). 

To  determine  if  the  continuous  poleward  flow  around  Point  Conception  is 
dependent  on  the  irregular  coastline  geometry  of  the  Southern  California  Bight,  a  second 
experiment  was  run.  The  second  experiment  is  identical  to  the  basic  experiment  except 
that  a  straight  coastline  south  of  Point  Conception  (~35°  N)  is  used  (see  Figure  18  for  the 
coastline  geometry).  By  comparing  and  contrasting  the  results  of  both  experiments,  we 
gain  a  better  understanding  of  the  importance  of  irregular  coastline  geometry  south  of 
Point  Conception  in  the  formation  of  surface  poleward  flow. 

As  expected,  the  initial  oceanic  response  to  the  combination  of  thermohaline 
gradients  and  wind  forcing  results  in  a  poleward  (equatorial)  flow  in  the  northern 
(southern)  portion  of  the  model  domain.-  By  spring  with  the  arrival  of  the  North  Pacific 
High  along  the  west  coast,  surface  equatorward  flow  extends  through  the  model  domain. 
A  coastal,  poleward  undercurrent  develops  below  the  surface  equatorward  current, 
initially  in  the  equatoward  end  of  the  model  domain  and,  during  the  upwelling  season, 
along  the  entire  coast.  During  the  summer  and  early  fall  as  the  coastal  jet  and 
undercurrent  become  fully  established,  the  currents  become  unstable  and  develop 
meanders,  eddies,  and  filaments  in  the  equatorward  end  of  the  model  domain  (e.g..  Figure 
18). 
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Similar  to  the  basic  experiment,  in  the  fall  and  winter,  the  coastal  poleward 
undercurrent  subsequently  strengthens  and  shoals  (e.g..  Figures  19a-19c).  By  September 
(Figure  19a),  poleward  flow  of  ~7.5  cm/s,  underlying  equatorward  flow  at  the  surface,  is 
generated  with  its  core  centered  below  ~500  m  depth.  By  October  (Figure  19b)  the  core 
of  the  undercurrent  strengthens  to  ~12.5  cm/s  below  ~5C)0m  depth.  By  November  (Figure 
19c),  a  coastal  poleward  current  of  ~7.5  cm/s  is  present  at  the  surface. 

Unlike  the  basic  experiment,  a  continuous  inshore  surface  poleward  flow  around 
Point  Conception  during  the  fall  and  winter  (or  any  other  season)  is  not  generated.  For 
example,  south  of  Point  Conception,  from  ~33.5°  N  to  ~34.5°  N,  inshore  surface 
poleward  flow  is  present  (e.g.,  Figure  20a).  Off  Point  Conception,  from  -34.6°  N  to 
-35.1°  N,  there  is  inshore  surface  equatorward  flow  (e.g..  Figure  20b).  North  of  Point 
Conception,  from  -35.2°  N  to  -36°  N  inshore  surface  poleward  flow  is  again  discernible 
(e.g..  Figure  19c). 

A  comparison  of  the  two  .  experiments  in  relation  to  surface  poleward  flow  (i.e., 
the  DC)  shows  that  both  model  simulations  generate  a  strengthening  and  shoaling  of  the 
poleward  undercurrent  nearshore,  while  only  the  basic  experiment  generates  a  continuous 
inshore  surface  poleward  flow  around  Point  Conception.  The  effect  of  irregular  coastline 
geometry,  specifically  the  Southern  California  Bight,  is  thus  shown  to  play  an  important 
role  in  generating  continuous  surface  poleward  flow  around  Point  Conception.  Based  on 
these  results,  both  the  shoaling  of  the  CUC  and  the  SCC  rounding  Point  Conception 
contribute  to  the  formation  of  the  DC. 
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Overall,  results  from  the  model  simulations  showed  that  the  combined  effects  of 


seasonal  wind  forcing  and  thermohaline  gradients  along  with  a  realistic  coastline  led  to 
the  generation  of  all  of  the  major  currents  of  the  CCS  as  well  as  the  formation  of 
filaments,  meanders,  and  eddies.  Tbe  eddies  were  shown  to  be  generated  from 
instabilities  of  the  CC  and  CUC  via  barotropic  and  baroclinic  instability  processes.  The 
relatively  large  magnitudes  of  MKE  compared  to  those  for  EKE  were  consistent  with  the 
hypothesis  that  eddies  are  driven  by  the  mean  flow  instead  of  the  mean  flow  driving  the 
eddies. 

With  the  larger  model  domain,  the  model  results  showed  that  south  of  Cape 
Blanco  the  equatorward  jet  was  a  continuous  feature,  meandering  offshore  and  onshore, 
and  divided  coastally-influenced  water  from  water  of  offshore  origin.  The  area  off 
southern  Baja  was  shown  to  be  a  highly  dynamic  environment  for  meanders,  filaments, 
and  eddies,  while  the  area  off  Point  Eugenia,  which  represents  the  largest  coastline 
perturbation  along  the  Baja  peninsula,  was  shown  to  be  a  persistent  cyclonic  eddy 
generation  region.  Lastly,  the  results  of  the  sensitivity  study  showed  that  both  the  SCC 
rounding  Point  Conception  and  the  shoaling  of  the  poleward  undercurrent  play  important 
roles  in  generating  the  DC  in  the  fall. 
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Model  Domain 


Longitude  (Degrees  West) 


Figure  1.  Model  domain,  coastline,  and  generalized  classical  circulation  of  the  California 
Current  System  (CCS).  The  domain  is  bounded  by  22.5°  N  to  47.5°  N,  107.5°  W  to 
132.5°  W.  The  broad,  slow  surface  equatorward  California  Current  (CC)  overlies  the 
narrow,  poleward  California  Undercurrent  (CUC).  Surface  poleward  flows  include  the 
Davidson  Current  (DC)  north  of  Point  Conception,  and  the  Southern  California  Eddy 
(SCE)  and  Southern  California  Countercurrent  (SCC)  south  of  Point  Conception. 
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Figure  3.  Time  series  plots  of  monthly  temperature  fields  used  as  seasonal  forcing  in  the 
basic  simulation.  The  symbol  represents  data  at  22.5°  N,  132.5°  W,  while  the  ‘+’ 
symbol  represents  data  at  47.5°  N,  132.5°  W  for  13  m  depth. 
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Figure  4.  Day  45  density  contours  and  velocity  vectors  at  13  m  depth  for  (a)  the  northern 
half  and  (b)  the  southern  half  of  the  model  domain.  Contour  interval  is  0.1  g/cm^; 
maximum  velocity  vector  is  50  cm/s. 
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Figure  5.  Day  195  density  contours  and  velocity  vectors  for  the  northern  (a)  and 
southern  (b)  portions  of  the  model  domain  at  13  m  depth.  Contour  interval  is  0.2  g/cm^; 
maximum  velocity  vector  is  100  cm/s. 
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Figure  7.  Cross-shore  sections  of  meridional  velocity  (v)  at  (a)  25°  N  (off  Cape  San 
Lazaro)  and  (b)  43°  N  (off  Cape  Blanco)  on  day  285.  The  contour  interval  is  5  cm/s  (10 
cm/s)  for  poleward  (equatorward)  flow  in  (a)  and  2.5  cm/s  (5  cm/s)  for  poleward 
(equatorward)  flow  in  (b). 
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Figure  8.  Day  285  velocity  vectors  for  the  northern  (a)  and  southern  (b)  portions  of  the 
model  domain  at  13  m  depth.  Maximum  velocity  vector  is  100  cm/s. 


45.0°N  H 


U  41.0"N  H 
Q  j 


3S.0®N 


37.0‘^N 


35.0°N 


If  ! 


I'' . -'®^-*- - 


icpeMendccIho 


2  W  '«^ 


’?0Rt  Wcf, 

“4^'; 

'Va'''v 


131.C°W  129.0®W 


i  I  i 

127.0°W  125,0‘^W 

LONGITUDE 


I  I  I 
123.0°W  121.C°W 


UUIMVjI  i  uuc. 

MEAN  KINETIC  ENERGY  ((cin/sec)**2) 


Figure  9.  Horizontal  maps  at  13  m  depth  of  mean  kinetic  energy  (MKE)  for  the  northern 
portion  of  the  domain  in  model  year  3  averaged  for  the  months  of  (a)  July  and  (b) 
October.  Contour  interval  is  200  cmVs^. 
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Figure  10.  Horizontal  maps  at  13  m  depth  of  eddy  kinetic  energy  (EKE)  for  the  northern 
portion  of  the  domain  in  model  year  3,  averaged  for  the  months  of  (a)  July,  (b) 
September,  and  (c)  October.  Contour  interval  is  25  cm^/s^. 
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Figure  11.  Horizontal  maps  at  13  m  depth  of  mean  kinetic  energy  (MKE)  for  the 
southern  portion  of  the  domain  in  model  year  3  averaged  for  the  months  of  (a)  July  and 
(b)  October.  Contour  interval  is  200  cm^/s^. 
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Figure  12.  Horizontal  maps  at  13  m  depth  of  eddy  kinetic  energy  (EKE)  for  the  southern 
portion  of  the  domain  in  model  year  3,  averaged  for  the  months  of  (a)  July  and  (b) 
October.  Contour  interval  is  25  cm^/s^. 
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Figure  13.  Cross-shore  sections  of  salinity  on  day  924  (mid-July)  of  model  year  3  at  26° 
N  (south  of  Point  Eugenia).  The  contour  interval  is  0. 1  psu. 


59 


I  lid  JO 


MEAN  MERID  VELOCITY  (cm/sec) 


Figure  14.  Cross-shore  sections  of  mean  v  at  25°  N  (off  Cape  San  Lazaro)  for  (a) 
January,  (b)  April,  (c)  July,  and  (d)  October  of  model  year  3.  The  contour  interval  is  2.5 
cm/s  (5  cm/s)  for  poleward  (equatorward)  flow  in  (a,  b,  d)  and  5  cm/s  for  both  in  (c). 
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Figure  15.  Cross-shore  sections  of  mean  v  for  January  of  model  year  3  at  (a)  36.5°  N  (off 
Point  Sur)  and  (b)  34.6°  N  (off  Point  Conception).  The  contour  interval  is  2.5  cm/s  (5 
cm/s)  for  poleward  (equatorward)  flow  in  (a)  and  4  cm/s  (5  cm/s)  for  poleward 
(equatorward)  flow  in  (b). 
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Figure  16.  Cross-shore  sections  of  mean  v  at  32°  N  for  (a)  September,  (b)  October,  and 
(c)  November  of  model  year  3.  The  contour  interval  is  2.5  cm/s  (5  cm/s)  for  poleward 
(equatorward)  flow. 
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Figwre  17.  Cross-shore  sections  of  mean  v  for  November  of  model  year  3  at  (a)  33.5°  N 
(south  of  Point  Conception),  (b)  34.6°  N  (off  Point  Conception),  and  (c)  35.5°  N  (north  of 
Point  Conception).  The  contour  interval  is  2.5  cm/s  (5  cm/s)  for  poleward  (equatorward) 
flow. 
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Figure  18.  The  straight  coastline  south  of  Point  Conception  used  in  experiment  2  is 
depicted  by  the  shaded  area.  Temperature  contours  and  velocity  vectors  in  experiment  2, 
for  day  285  (mid-October).  The  contour  interval  is  1°  C;  maximum  velocity  vector  is  100 
cm/s. 
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Figure  19.  Cross-shore  sections  of  mean  v  of  e> 
September,  (b)  mid-October  and  (c)  mid-Novembei 
cm/s)  for  poleward  (equatorward)  flow. 


DEPTH  :  13m 


45.0°  N 


15.4 


15.2 


14.8 


14.6 


14.4 


14.2 


35.0°N 


13.8 


13.6 


13.4 


127.0°W  125.0°W  123.0'’W  121.0°W 

LONGITUDE 


131.0°W  129.0°W 

KINETIC  ENERGY  TRANSFER  ((cm/sec)**2) 


Plate  1.  Horizontal  maps  at  13  m  depth  of  mean  kinetic  energy  transfers  (a,  c)  and 
potential  energy  transfers  (b,  d),  averaged  over  days  180-195  in  the  northern  portion  of 
the  domain  (a,  b)  and  the  southern  portion  of  the  domain  (c,  d).  Contour  interval  is  0.2 
cm^/s^  for  (a),  2.0  cmVs^  for  (b,  c),  and  5.0  cm^/s^  for  (d). 
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Plate  2.  Horizontal  maps  at  13  m  depth  of  mean  kinetic  energy  transfers  (a,  c)  and 
potential  energy  transfers  (b,  d),  averaged  over  days  270-285  in  the  northern  portion  of 
the  domain  (a,  b)  and  the  southern  portion  of  the  domain  (c,  d).  Contour  interval  is  5.0 
cm^/s^  for  (a),  (c),  and  (d),  and  2.0  cmVs^  for  (b). 
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Plate  3.  Mean  temperature  and  velocity  vectors  at  13  m  depth  for  (a)  January,  (b)  April, 
(c)  July,  and  (d)  October  of  model  year  3.  Contour  interval  is  0.1°  C;  maximum  velocity 
vector  is  100  cm/s. 
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Table  1.  Values  of  Constants  Used  in  the  Model 


Constant 

Value 

Definition 

To 

278.2®K 

Constant  Reference  Temperature 

^0 

34J 

Constant  Reference  Salinity 

Po 

Density  of  Sea  Water  At  Tq  and  Sq 

a 

2.4  X  IQ"*  (°K)-’ 

Thermal  Expansion  Coefficient 

7.5  X  10-^ 

Saline  Expansion  Coefficient 

K 

10 

Number  of  Levels  In  Vertical 

Ax 

Cross-Shore  Grid  Spacing 

Ay 

Alongshore  Grid  Spacing 

H 

Total  Ocean  Depth 

At 

800  s 

Time  Step 

fo 

Mean  Coriolis  Parameter 

g 

Acceleration  of  Gravity 

2  X  lO^^cm^s’^ 

Biharmonic  Momentum  Diffusion  Coefficient 

Ah 

2  xlO’^cmS'* 

Biharmonic  Heat  Diffusion  Coefficient 

Km 

0.5  cm^  s'^ 

Vertical  Eddy  Viscosity 

Kh 

0.5  cm^  s'^ 

Vertical  Eddy  Conductivity 
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